I.  INTRODUCTION 


With  thfi  increasing  in'.portancc  in  modern  vmrfare  of 
infrared  laser  weaponry,  guidance  systems,  and  thermal  imaging 
devices,  it  has  become  more  and  more  important  that  v/e  under¬ 
stand  techniques  for  obscuration  by  particulate  systems,  i.e. 
smoke  screens.  A  light  beam  is  attenuated  by  the  two  processes, 
absorption  and  scattering,  t}ie  sum  of  those  being  called  extinc¬ 
tion.  A  desirable  property  of  the  screening  smoke  is  high 
extinction  efficiency  per  unit  volume  or  per  unit  mass  of  the 
material  (depending  on  the  exact  application) .  In  order  to 
attenuate  CO2  laser  radiation  one  requires  a  particulate  system 
having  large  extinction  coinciding  with  the  CO2  laser  emission 
band  wavelengths  of  approximately  9.6  and/or  10.6  )jm.  Band 
widths  need  not  be  large  in  this  case.  For  the  purpose  of 
screening  against  therjaal  im^fing  devices,  which  may  utilize  the 
entire  atmospheric  window  region  betv;een  about  8  and  14  pm 
v;avolength,  a  considerably  wider  particulate  extinction  band 
is  called  for.  One  mig)it  at  first  think  that  the  best  approach 
would  be  to  use  an  aerosol  v;hich  would  give  higli  extinction 
throughout  the  infrared.  It  is  possible,  hov;evcr,  to  achieve 
much  higher  extinction  in  a  narroT^band  than  in  a  wide  band, 
as  will  be  shown  in  this  report. 

A.  Purpose  of  this  work 

In  order  to  successfully  tailor  the  aerosol  to  the  specific 
need  we  should  be  able  to  engineer  the  optical  properties  of 
the  particulates,  that  is  to  calculate  the  performance  of 
various  systems  and  to  produce  the  one  that  most  closely  fits 
the  desired  characteristics.  In  order  to  accomplish  this 
optical  engineering  of  aerosols  the  following  two  things  are 
necessary : 

(1)  A  catalogue  of  the  wavelength-dependent  optical  v 
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'  constants  for  various  solirls  and  liquids. 

(2)  Adequate  theoretical  representation  of  hov/  extinction 
is  related  to  such  factors  as  particle  size,  shape,  and  coating. 


Groat  difficulties  have  boon  present  in  both  of  these  areas. 
Regarding  the  all  important  optical  constants,  there  is  no 
extensive  compilation  available.  This  is  true  in  spite  of  the 
fact  that  there  is  a  large  amount  of  literature  on  certain 
optical  characteristics  such  ns  transmission.  Optical  constants, 
hov;evor,  are  usually  determined  only  for  liomogcneous  samples 
such  as  single _£rystals  and  bul}:.O.icjuids .  Unfortunately,  most 
solids  have  not  been  produced  as  single  crystals,  but  are 
commonly  available  only  in  povnlcxcd  form.  A  technique  to 
permit  accurate  determination  of  optical  constants  for  pov/dered 
materials  would  provide  a  much  needed  means  for  filling  this 
gap  in  our  basic  physical  data _bank.  The  main  goal  of  this 
work  has  been  to  attempt  development  of  such  a  technique. 

In  the  course  of  the  technique  development  wo  have  had 
to  deal  critically  with  point  (2)  above,  the  theory  of  sjnall 
particle  extinction.  Some  progress  in  the  theoretical  treatment 
of  shape  effects  has  been  achieved  as  a  necessary  part  of  the 
optical  constants  technique  development.  The  purpose  of  this 
report  is  to  describe  the  development  of  the  technique  for 
determining  optical  constants  of  powdered  solids  and  to  summa¬ 
rize  the  procedure.  Application  of  the  method  to  four  solids 
with  well-known  optical  constants  provides  a  means  for  eval¬ 
uating  the  success  of  the  procedure.  As  a  by-product,  some 
insights  we  have  gained  into  the  larger  question  of  how  to 
produce  maximum  extinction  will  be  discussed. 


B.  Organization  of  this  report 

Organization  of  this  report  will  bo  in  the  following 
way.  In  part  II  we  discuss  methods  for  determining  optical 


constants  starting  with  the  well-tested  methods  used  on  smooth, 
bulk  materials.  This  enables  us  to  present  a  yiev/point  on 
'how  to  treat  the  similar  but  much  more  difficult  problem  in 
powdered  samples.  Difficulties  in  previous  attempts  at 
determining  optical  constants  for  powders  are  discussed,  and  the 
philosophy  of  the  present  v;ork  is  presented,  emphasizing  the 
special  need  for  better  understanding  of  absorption  by  irregu¬ 
larly  shaped  particles.  In  part  III  we  point  out  why  it  is 
desirable  to  limit  both  experiments  and  theory  to  the  Rayleigh 
regime.  Part  IV  develops  a  theory  for  treating  irregular 
shapes  which,  v/hile  not  perfect,  gives  a  significant  im¬ 
provement  over  the  theory  for  spheres.  Part  V  summarizes 
the  experimental  techniques  for  measuring  extinction  by 
small  particles  and  compares  the  measurements  with  calcula¬ 
tions  biiscd  on  the  shape  distribution  theory  of  IV  for  the 
cases  of  four  solids  with  well-knov.’n  optical  properties. 

Part  VI  presents  our  attempts  to  invert  the  small  particle 
extinction  data  for  the  four  "knov;n"  solids  and  discusses 
the  degree  of  success.  The  larger  question  of  how  does  one 
produce  maximum  extinction  by  moans  of  absoriDtion  by  particles 
in  the  10  yrn  spectral  region  is  retairned  to  in  part  VJT 
Several  insights  into  this  question  gleaned  from  the  present 
combination  of  theory  and  experiment  are  coinmentcd  on.  Finally, 
part  VIII  summarizes  the  conclusions  of  this  work  and  makes 
suggestions  for  further  work. 

II.  DISCUSSION  OF  SOME  METHODS  FOR  DETERMINING 

OPTICAL  CONSTAr.'TS 

The  basic  optical  properties  of  matter,  such  as  absorp¬ 
tion,  reflection,  and  scattering.,  are  determined  by  the 
wavelength-dependent  complex  optical  constants 

m{X)  -  n  +  ik 

where  n  is  the  index  of  refraction  and  k  is  the  extinction 


coefficient  (although  this  terminology  will  be  avoided  in 
this  report  because  of  possible  confusion  with  sma 1 1  particle 
extinction) .  There  are  numerous  common  misconceptions  con¬ 
cerning  optical  constants,  such  as  that  n  must  be  greater  than 
one  to  avoid  violation  of  Einsteinian  relativity  —  not  true; 
such  as  that  n  may  he  considered  as  approximately  constant 
as  k  varies  not  alv;nys  true.  Those  misconceptions,  both 
of  which  play  a  part  in  this  discussion,  probcibly  arise 
because  toclmical  people  have  been  trained  by  their  elementary 
optics  teaching  that  n  is  essentially  constant  (  -1.5)  with 
only  a  small  dispersion  tov;ard  the  blue  that  increases  n  by 
a  fev;  percent,  while  J:  is  usually  zero  but  may  vary  a  little 
in  colored  materials.  These  statements  are,  after  all,  true 
for  glasses  in  the  visible,  which  is  what  traditional  optics 
is  all  al:iout .  The  truth  is,  hov;evcr,  that  in  the  more 
remote  v;avolength  regions  of  the  ultraviolet  and  the  infrared 
all  solids  have  regions  in  which  both  n  and  k  vary  over 
v/ide  ranges.  Very  commonly  n  is  less  than  1,  k  is  greater 
than  n,  and  both  vary  over  orders  of  magnitude  v/ithin  a 
narrow  v/ave length  region.  Examples  of  this  behavior  will 
come  up  in  this  report. 

‘  In  addition  to  the  indices  of  refraction,  optical 
properties  can  be  summarized  in  terms  of  complex  dielectric 
functions 

e(X)  =  e  +  ie, 

1  *  • 

The  complex  indices  of  refraction  and  the  complex  dielectric 
functions  are  directly  related  in  the  following  way; 


G ,  =  n*  -  k* 

Gg  =  2nk 


(1) 


Hence  cither  pair  of  quantities  can  bo  said  to  completely 
describe  the  optical  properties  of  a  substance,  and  we  shall 
frequently  use  the  term  optical  constants  to  refer  to 
either  set. 
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A.  Bulk  optical  constants 

The  most  direct  determination  of  optical  constants 
is  accomplished  by  measurement  of  transmission  through  a 
slab  of  the  material  and  reflectance  from  its  surface. 

For  v/avolength  regions  and  for  solids  of  concern  in  this 
v;ork,  however,  <'())Sorption  in  so  great  that  bulk  solids  cannot 
be  produced  in  sufficiently  thin  slices  that  any  measured 
light  can  penetrate  through  the  slab.  One  must  necessarily 
resort  to  pure  reflectance  tec])niques.  Since,  however,  one 
measurement  cannot  determine  the  two  desired  quantities,  the 
usual  procedure  is  to  measure  reflcct^lnce  under  two  different 
conditions  (tv.’o  polarisation  states  at  non-normal  incidence 
or  two  different  angles)  or  to  measure  reflectance  over  an 
extended  v/avelength  range  coupled  with  an  analysis  that  proper¬ 
ly  relates  n  and  k.  It  is  this  last  technique  that  v;e  shall 
present. 

The  normal  incidence  reflectance  R  of  a  surface  of 
material  charactcrix.ed  by  opticiil  constants  n  and  k  is  given 
by  the  Fresnel  equation 


R  = 


(n-1)^  +  k^ 
(n+1)^  + 


(2) 


v;here  the  material  is  assumed  to  be  in  air  with  index  of 
refraction  taken  to  be  1.  When  tlie  reflectance  of  certain 
solids  is  measured  in  the  infrared  there  is  a  region  of 
wavelengths  for  vdiich  R  is  very  close  to  1.  An  example  is 
shown  in  the  a  part  of  figure  (1)  for  the  solid  silicon  carbide. 
The  historical  terminology  for  the  higlily  reflecting  region  is 
"reststrahlcn''  region,  since  multiple  reflection  of  polychro¬ 
matic  light  from  such  a  solid  results  in  ’'residual  radiation" 
at  these  wavelengtlis .  The  reststrahlcn  region  is  caused  by 
strong  interaction  of  light  with  the  ionic  vibration  modes  of 
the  crystal.  Almost  all  solids  and  liquids  have  vibrational 
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Figure  1:  (a)  Measured  reflectance,  (b)  derived  optical 

constants,  and  (c)  derived  dielectric  functions  for  a-SiC 
from  Spitzer,  Kleinman,  cuid  Walsh  (1959).  The  circles  in 
the  reflectance  graph  represent  some  of  the  data  points  and 
the  continuous  curves  in  all  three  graphs  are  calculated  from 
the  one  oscillator  fit  to  the  data  with 

Ej"  6.7  w,  -  793  (i)p  ■  1443  y  »  4.76  . 


absorption  bands  in  the  infrared,  but  the  strengths  of  the  bands 
may  not  be  groat  enough  to  produce  a  reflectance  approaching  1 
as  in  the  roststrahlen  case.  It  is,  hov/over,  the  strong 
vibrational  interaction  that  ultimately  gives  rise  to  the 
strong  small  particle  absorption  bands  that  are  of  chief  concern 
in  this  work.  Hence  our  attention  is  focused  primarily  on  bands 
such  as  those  of  figure  (la) .  For  such  cases  it  has  been  found 
that  a  charged  liarmonic  oscillator  model  of  the  crystal  vibra¬ 
tions  and  their  interaction  with  the  electromagnetic  field  gives 
a  quite  accurate  description  of  the  optical  effects  in  the 
bulk  solid.  For  a  solid  composed  of  N  identical  oscillators 
per  unit  volume  the  complex  dielectric  function  is  (Hodgson 
1970,  for  example) 


where  =  Ho^/m  i.s  the  plasma  frequency  of  the  ions  of  mass  m. 
(Oj  is  the  vibrational  frc<rucncy  related  to  the  effective  spring 
constant  and  the  mass,  and  y  is  the  damping  factor.  is 

the  dielectric  function  at  high  frequencies  compared  to  t)ie 
oscillator  frequency.  The  high  frequency  dielectric  constant 
is  normally  due  to  the  effect  of  electronic  oscillators  in 
the  ultraviolet.  If  there  arc  several  types  of  oscillators, 
due  for  example  to  different  bonds  bctv;een  ions,  equation  (3) 
can  be  written  as  a  sum  over  the  several  types  j , 


The  multiple  oscillator  model  forms  the  basis  of  a  common 
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technique  for  analyzing  normal  incidence  reflectance  data  for 
bulk  optical  constants  of  the  solid.  The  results  are  illus¬ 
trated  in  figure  (])  for  the  solid  silicon  carbide.  Values  of 
the  paranotors  c  ^  ty ,  Wj,  and  (jj^are  assumed  (guessed)  and  used  in 
equation  (3)  to  calculate  c,(X)  and  C2(X)  ,  which  are  converted 

to  n(X)  cxnd  k(X)  using  equations  (1).  Indices  of  refraction 
are  then  used  in  equation  (2)  to  calculate  P.{X)  .  Q'he  parameters 
arc  varied  until  the  calculated  reflectance  gives  a  good 
match  to  the  laeasurcd  reflectance  spectrum.  Once  the  fit  is 
achieved  the  oscillator  parameters  then  give  ej(X)  and  E2(X) 
or  n(X)  and  k(X)  as  desired.  The  full  curves  in  figure  (la) 
are  the  calculated  values  an  compared  v;ith  measured  reflectance 
values  denoted  by  circles.  Optical  constants  in  (l^-*)  and  (Ic) 
are  the  inferred  values  based  on  the  best  fit. 

If  the  solid  is  anisotropic  the  technique  can  still  bo 
employed  provided  the  sample  is  cut  and  oriented  properly  and 
polarized  liglit  is  used.  A  classic  example  of  optical  constants 
determination  for  an  anisotropic  solid  using  the  multiple 
oscillator  expression  (i)  is  the  v.’ork  on  crystalline  quartz 
of  Spitzer  and  Kleinman  (1961) .  Figure  (2)  shows  the  measured 
reflectance  for  the  tv;o  rnajor  crystalline  directions  as  solid 
lines  ^lnd  the  fits  to  the  data  achieved  using  the  multiple 
oscillator  theory.  Their  derived  parameters  and  the  optical 
constants  resulting  from  them  v/ill  by  used  in  ovr  work  discussed 
later  in  this  report. 

There  are  several  reasons  for  our  presenting  tliese  details 
of  the  technique  for  determining  optical  constants  by  using 
the  oscillator  model.  Wo  v;ill  use  an  analogous  method  in  our 
task  of  determining  n  and  k  for  pov;dorcd  materials.  Secondly, 
we  wish  to  assure  the  reader  that  the  oscillator  model, 
although  very  simple,  has  convincing  validity  V7hcn  applied  to 
solids  such  as  v/c  are  most  interested  in.  Finally,  this 
commonly  used  method  for  "measuring"  optical  constants  em¬ 
phasizes  the  necessity  of  interposing  theory  (in  this  case 
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Figure  2  Reflectivity  of  oriented  single  crystal  quartz 
for  the  two  major  polarization  directions.  Data  are 
shown  by  the  points  and  the  theoretical  fit  by  the 
solid  lines. 


Lorentz  oscillator  theory  plus  Fresnel's  equation)  between  the 
measurements  and  the  derived  quantities.  Diagramatically , 


Measurements  Theory 


Another  popular  technique  of  analyzing  normal  incidence  reflec¬ 
tance  data  for  extracting  optical  constants  is  by  means  of 
so-called  Kramers-Kronig  relations.  This  technique  has  been 
discussed  in  many  places  (see  for  example  Stern  1964).  The 
diagram  is 

Measurements  Theory 


2 

It  is  frequently  true,  as  in  the  above  situations,  that  the 
quantity  "measured"  is  actually  an  inferred  quantity  based  on 
the  actual  measurement  of  a  related  quantity  coupled  with 
some  sort  of  connecting  link,  the  tlieory.  Assuming  that  the 
measurements  are  perfectly  correct,  the  resxilts  can  still  be 
accurate  only  if  the  intermediary  theory  is  accurate  and  applies 
accurately  to  the  experimental  situation.  In  the  case  of  the 
above  mentioned  determinations  of  optical  constants  for  quartz 
and  for  SiC,  we  have  great  confidence  in  Fresnel's  equations, 
the  Kramers-Kronig  relations,  and  the  Lorentz  oscillator 
equations.  Furthermore,  single  crystals  such  as  were  used  in 
those  examples  can  be  cleaved  or  polisliod  to  insure  that  the 


normal  incidence  reflection  equation  is  accurately  obeyed  in 
the  experimental  measurements. 

A  more  subtle  feature  of  the  two  bulk  optical  constants 
methods  is  that  n  and  k  are  correctly  related  to  one  another. 

It  is  an  important  point  in  this  work  that  n  and  k  (or  Sj  and  c^) 
are  not  independent  quantities  but  are  real  and  imaginary 
parts  of  the  same  complex  quantities.  As  such  they  are 
inseparably  related  to  one  another.  The  relationship  of  Sj  and  Cj 
to  one  another  can  be  expressed  in  fact  by  another  one  of 
the  Kreimers-Kronig  relations. 

ej({o)  =  1  +  -  I  -  du' 

j  (o'  -  (0^ 

(5) 

(oe  j  (co ' ) 

^2^^)  -  i-  [ - 

^  TT  J  •  2 

'  (0  -  (0 


Any  method  that  determines  a  wavelength-dependent  e j (or  n) 
by  one  technique  and  a  wavelength  dependent  Cj (or  k)  by  an 
unrelated  technique  that  does  not  properly  aoknov/ledge  the 
inter-relationship  of  the  quantities  should  bo  used  with 
caution. 

B.  Small  particle  optical  constants 

Several  attempts  have  been  made  to  determing  optical 
constants  of  powdered  or  otherwise  inhomogeneous  materials. 
Volz  (1972),  for  example,  has  presented  "measured"  optical 
constants  for  atmospheric  dust  particles,  Perry  et.  al.  (1972) 
have  given  optical  constants  for  lunar  soil  samples,  and 
Penman  (1976a, b)  has  inferred  optical  constants  from  minerals 
and  meteoritic  material  bearing  a  possible  relation  to  inter¬ 
stellar  grains.  The  subject  of  our  present  study  is  (fortu¬ 
nately)  a  little  simpler  than  most  of  these  tasks.  We  are 
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dealing  with  only  one  chemical  component,  such  as  SiC  or  CaCO^, 
in  pov/dered  form,  as  opposed  to  a  rock,  a  lunar  soil  sample,  a 
collection  of  atmospheric  dust,  or  a  meteorite,  each  of  which 
are  multicomponent  mixtures.  Some  of  the  work  by  Volz  and 
some  by  Penman  v/cre  concerned  v/ith  single  component  powders, 
hov/ever,  and  will  be  discussed  briefly.  A  recent  paper  by 
Jennings  et.  al.  (1979)  has  been  critical  of  the  techniques  for 
determining  optical  constants  of  inhomogeneous  materials, 
particularly  for  atmospheric  particulates.  We  now  discuss  .some 
of  the  difficulties  in  relation  to  previous  methods  of  deter¬ 
mining  optical  constants  and  to  our  ovm  proposed  method.  Our 
conclusions  are  that,  despite  difficulties,  there  is  hope  of 
extracting  reasonably  reliable  optical  constants  from  one- 
component,  powdered  materials.  The  discussion  of  difficulties 
emphasizes  experimental  and  theoretical  problem.s  to  which 
particular  attention  must  bo  paid. 

The  metliod  v/hich  Volz  (1972)  has  used  for  atmospheric 
particulates  is  summarized  in  our  diagrammatic  way  as  follows; 


Measurements 


'Dieory 


Transmission 

Bulk 

and  Reflectance 

transmission 

vs  X 

for  powder  in  KBr 

& 

ref lectance 

and  pure  pov;der 

theory 

Volz  mixed  a  .small  quantity  of  sample  pov,»der  into  KBr,  formed  a 
pellet  and  measured  tran.smission.  The  transmission  (T)  was 
analyzed  using  the  equation 

T  =  exp(-  --^t)  (6) 


where  X  is  the  wavelength  and  t  is  the  effective  thickness  of 
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the  sample  pov/der.  This  equation  treats  the  sample  powder  as 
a  homogeneous  slab.  The  procedure  has  been  justly  criticized 
by  Jennings  et.  al.  for  not  taking  into  account  certain  small 
particle  properties  such  as  scattering  and  the  size  dependence 
of  extinction.  Volz  continues  his  determination  of  n  by  using 
the  k's  derived  from  equation  (6)  along  with  measured  normal 
incidence  reflectance  R  of  compressed  powder  (with  no  KBr)  in 
equation  (2) .  In  many  cases  powders  cannot  be  compressed  to 
anything  like  bull:  solid  surfaces  so  that  the  bulk  reflectance 
equation  is  not  applicable.  Even  if  (6)  v;ore  perfectly  appli¬ 
cable  and  if  there  were  no  scattering  effects  and  size  depen¬ 
dence  of  extinction,  the  serious  errors  introduced  by  treating 
extinction  as  a  bulk  property  in  determining  k  would  carry  over 
to  introduce  serious  errors  in  n. 

We  have  tried  Volz's  method  in  an  attempt  to  derive 
optical  constants  of  MgO.  Thcr  results  were  so  far  off  that 
no  purpose  seemed  to  be  served  by  shov.'ing  them.  The  reason  for 
the  severe  difficulty  in  this  technique  as  applied  to  solids 
with  strong  infrared  bands  is  as  follows.  Small  particles  of 
a  solid  have  resonant  absorption  modes  at  different  wavelength 
than  does  the  bulk  solid,  but  this  does  not  im.ply  that  the 
optical  constants  of  the  small  particle  and  the  bulk  are  differ¬ 
ent.  The  shift  in  absorption  bands  is  perfectly  calculable  from 
bulk  optical  constants.  For  spheres  the  small  particle  mode 
occurs  at  the  v;avelength  for  which  e,  =  -2c  e„  is  the  dielec- 
trie  function  for  the  surrovmding  medium,  in  this  case  Klir . 

The  bulk  solid  has  its  resonant  mode  at  the  frequency  where 
Cj  is  a  maximum.  Thus  if  equation  (6)  is  employed,  treating 
the  small  particles  as  bulk,  the  inferred  c^  (or  k)  is  sub¬ 
stantially  displaced  from  its  true  position.  Introduction  of 
this  highly  erroneous  k  into  the  reflectance  equation  then  gives 
rise  to  even  worse  values  of  n.  Without  desiring  to  be  critical 
we  have  mentioned  this  technique  to  emphasize  the  errors  that 
can  accrue  V7hcn  proper  account  is  not  taken  of  small  particle 
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as  opposed  to  bulk  properties. 

Penman's  (197G)  approach  has  been  to  measure  specular 
reflectance  at  almost-normal-incidence  for  his  heterogeneous 
solids  and  to  analyze  by  means  of  the  Fresnel  reflectance 
equation  (2)  and  Kramers-Kronig  analysis,  as  though  the  sample 
were  a  homogeneous  solid  with  a  smooth,  plane  surface. 


Measurements  Theory 


Normal  incidence 

Fresnel  refl. 

reflectance 
vs.  X 

+ 

Kramers-Kronig 

analysis 

This  is  the  exact  technique  that  has  been  used  often  for  bulk 
solids,  but  the  theory  is  not  strictly  applicable  to  the 
inhomogeneous  samples,  and  the  effect  of  optical  anisotropy 
is  not  taken  into  account. 

The  approach  we  outlined  in  the  proposal  for  this  con¬ 
tract  was  to  use  transmission  measurements  of  the  particulates 
as  a  dilute  suspension  in  a  KBr  pellet,  similar  to  what  Volz 
did.  We  proposed  to  analyze  the  transmission  measurements  by 
using  small  particle  absorption  and  scattering  tlieory  coupled 
with  a  multiple  oscillator  theory. 


Measurements 


Theory 


n,k 


The  merit  of  this  approach  as  contrasted  with  tJie  Volz  and 
Penman  approaches  is  that  it  takes  into  account  specifically 
the  small  particle  nature  of  the  sample.  Optical  constants 
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'  are  also  related  correctly  to  one  another  in  a  self  consistent 
way  by  the  oscillator  model.  Unfortunately,  our  early  work 
under  this  contract  proved  that  a  very  important  factor  is  not 
included  in  this  approach  —  the  shape  distribution  of  the 
particles . 

Although  it  is  commonly  believed  tliat  shape  effects  are 
not  important  in  the  Rayleigh  region,  especially  for  absorption, 
this  is  not  at  all  true  in  the  vicinity  of  very  strong  infrared 
resonances.  A  previev;  of  figure  (7)  ,  v;hich  v/ill  be  discussed 
in  more  detail  later,  will  shov;  the  reader  that  the  sphere 
theory  using  accurately  measured  bulk  optical  constants  of 
quartz  fails  badly  in  calculating  the  extinction  band  properties 
of  a  real  pov.'der  of  quartz  particles.  The  discrepancy  lies  in 
the  fact  that  absorption  modes  in  individual  particles  are 
quite  shape  dependent.  Because  of  th<-*  inadequacy  of  sphere 
theory  in  regard  to  shape  effects  v;e  have  made  considerable 
effort  in  attempting  to  develop  a  theory  to  correctly  treat 
irregular  shapes.  This  extensive  digression  was  not  planned 
at  the  beginning  of  the  vJor3:  but  turned  out  to  be  a  necessary 
deviation  from  the  original  plan.  Having  developed  a  more 
successful  small  particle  theory  for  treating  shape  effects, 
discussed  in  part  III,  our  method  for  inferring  optical  constants 
now  has  become 

Measurements  Theory 

n,k 


We  have  given  the  preceeding  discussion  in  order  to 
emphasize  that  one  does  not  measure  optical  constants  direct¬ 
ly;  one  performs  related  measurements  that  are  analyzed  with 


Transmission  of 
powder 
in  KBr 


Oscillator 

model 

+ 

shape  distribution 
tlieory 
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the  help  of  appropriate  theory  to  infer  n  and  k.  The  weakest 
part  of  the  process  for  small  particles  seems  to  be  the  theory 
box.  It  is  difficult  to  make  the  experiment  exactly  conform 
to  the  theory  (perfectly  isolated  spherical  particles,  for 
example) ,  and  it  is  difficult  to  construct  a  theory  that 
correctly  interprets  the  actual  experiment  (somewhat  clumpy 
collections  of  higlily  irregular  particles)  .  The  game  wo  play 
is  a  difficult  one,  and  the  results  can  hardly  be  expected  to 
be  as  accurate  as  single  crystal  results. 

III.  LIMITATION  OF  TilEORY  AND  EXPERIMI^NT  TO 
THE  RAYLEIGH  REGION 

Although  the  difficulties  in  inverting  experimental 
measurements  on  pov;ders  to  obtain  optical  constants  are  severe, 
great  simplifications  occur  if  only  particles  are  considered 
that  are  small  enough  for  the  Rayleigh  approximation  to  bo 
valid.  The  conditions  for  applici:vbility  of  the  Rayleigh  approx¬ 
imation  arc  that  the  particle  dimensions  (1)  are  small  com¬ 
pared  to  tine  v/avelcngth  (X)  and  that  the  phase  shift,  27rml/X, 
across  the  particle  is  small.  The  important  simplifications 
resulting  from  use  of  the  Rayleigh  theory  are  the  follov/ing: 

(1)  there  are  no  size  effects  in  volxime-normalized 
extinction . 

(2)  shape  effects,  which  arc  substantial  in  work  of  this 
sort,  can  be  simply  treated  thcoreticiilly . 

In  this  section  we  will  demonstrate  with  sphere  calculations 
how  the  size  dependence  goes  av;ay  for  small  radii,  then  show 
how  inadequate  is  the  spherical  tlicory  even  in  the  small 
particle  limit.  Finally,  v/e  preview  the  results  from  section 
IV  by  showing  the  improvement  a  Rayleigh  theory  of  shape  dis¬ 
tribution  can  provide  over  the  spherical  theory. 
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.  A.  Size  dependence 

The  only  shapes  of  arbitrarily  large  particles  for  which 
an  exact  solution  to  the  electromagnetic  boundary  value  problem 
is  possible  are  spheres  (the  Mie  theory)  and  infinitely  long 
cylinders.  By  using  Mie  theory  (see  van  de  Hulst  1957,  for 
example)  it  is  possible  to  determine  how  the  extinction  coeffi¬ 
cient,  normalized  per  unit  volume  of  solid,  varies  as  the  parti¬ 
cle  size  (radius)  changes.  Figure  (3)  shov7s  results  (Day  et.  al. 
1974)  of  Mie  calculations  based  on  measured  optical  constants 
at  infrared  v/avelengths  where  maximum  absorption  occurs  in  the 
vicinity  of  10  Urn.  Three  silicate  materials  were  chosen  for 
which  the  absorption  strengths  are  quite  large.  One  sees  that 
at  small  radii  (less  than  0.5  jjm)  the  volume  normalized  extinc¬ 
tion  is  constant.  For  radii  greater  than  about  1  ym  the  extinc¬ 
tion  is  extremely  size  dependent,  falling  quickly  as  radius 
increases.  This  decrease  of  extinction  with  radius  is  a 
saturation  phenomenon  which  occurs  when  absorption  is  so  effec¬ 
tive  that  the  inner  parts  of  tlie  sphere  are  hardly  influenced 
by  the  electromagnetic  fields.  As  the  radius  increases,  the 
volume  increases,  but  more  and  more  of  the  inner  volume  is 
ineffective;  thus  the  absorption  per  unit  volume  decreases. 

For  particles  in  the  1  to  10  ym  size  range  it  should  be  clear 
that  the  exact  size  distribution  must  be  known  and  taken  into 
account  in  order  to  infer  optical  constants  from  measured 
extinction.  For  particles  smaller  than  about  0.8  ym,  hov/ever, 
figure  (3)  shows  that  the  size  distribution  is  immaterial  since 
extinction  per  unit  volume  has  become  independent  of  size.  This 
is  the  consideration  demanding  that  we  experimentally  select 
out  the  sub-micron  size  fraction,  if  cit  all  possible,  for  doing 
our  extinction  measurements  on.  Jennings  et.  al.  (1979)  have 
called  attention  to  the  size  dependence  problem  in  analyzing 
optical  measurements  of  pov;dors  and  have  pointed  out  the 
necessity  of  knowing  the  size  distribution.  We  circumvent  the 
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CALCULATED  MASS  EXTINCTION  EFFICENCIES  ( 
AT  THE  IR  PEAK 


Figure  3 

Peak  mass  extinction  coefficient  (  »)  vs.  particle  radius  (a) 
for  amorphous  quartz,  as  well  as  crystalline  quartz  and  almandlne 
normalized  to  unity  for  spheres  with  radii  equal  to  O.lu. 

(  from  Day,  et.  al,,  1<J74)  ~ 
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»  need  for  knov/ing  the  size  distribution  by  selecting  only  par¬ 
ticles  small  enough  for  the  size  distribution  to  be  immaterial. 

In  order  to  confirm  the  size  dependence  and  to  develop 
techniques  for  selecting  a  suitably  small  size  fraction  we 
have  segregated  a  pov;dcr  of  quartz  into  several  size  fractions 
by  differential  settling  in  water.  Stokes’  law  for  the  rate 
of  settling  of  spherical  particles  of  radius  r  with  density 
falling  in  a  medium  with  density  and  viscosity  q  is 

(7)  " 

where  g  is  the  acceleration  due  to  gravity.  For  quartz  falling 
in  water  the  velocity  is  given  by 

V  =  3.6  X  10~^  r* 

where  v  is  in  cm/scc  and  r  is  in  ym.  Therefore  the  radius  of 
a  particle  can  be  determined  if  the  distance  traveled  and  the 
time  are  measured.  Experimentally,  approximately  10  grams  of 
a  finely  ground  quartz  powder  v;ere  mixed  v;ith  5  -  10  ml  of 
distilled  water.  A  100  ml  graduated  cylinder  whose  height  was 
marked  in  cm's  was  filled  nearly  to  the  top  with  distilled 
water.  The  dense  quartz  suspension  was  poured  into  the  cylinder 
and  allov;cd  to  settle  for  a  predetermined  length  of  time, 
after  which  successive  portions  of  the  fluid  v;ero  drawn  off 
from  the  top  and  placed  in  petri  dishes.  The  water  was 
then  allowed  to  evaporate,  leaving  behind  "sized"  particles. 

Extinction  measurements  were  made  for  a  few  different 
sized  particle  fractions  of  the  quartz  encapsulated  in  KBr 
pellets.  (See  section  V  for  a  more  detailed  description  of 
the  sample  preparation  and  measurement  techniques.)  Trans¬ 
mission  curves  for  three  of  the  size  fractions  in  KBr  are 
shown  in  figures  4,5,  and  6.  Figure  4  for  the  large  size 
fraction  also  shows  a  light  micrograph  of  the  particles,  con- 


Extinction  vs.  frequency  for  fused  silica  particles  with  ri2^ 
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f inning  that  they  are  inrlood  quite  irregularly  shaped.  Com¬ 
parison  of  the  three  spectra  including  the  derived  mass-nor¬ 
malized  extinction  coefficients  at  the  peak  shows  that  the 
shape  of  the  extinction  spectrum  for  largo  particles  is  quite 
different  from  the  spectrum  of  the  smaller  particles.  The 
peak  mass  extinction  coefficient  decreases  markedly  in  going  to 
the  smaller  sizes,  in  agreement  v/ith  the  discussion  of  figure 
(3)  .  Table  I  summarizes  the  i:)ea]:  mass  extinction  coefficients 
for  the  several  size  fractions  of  quartz. 


TABLE  I 

Mass  normalized  extinction  coefficients  for  quartz 
particles  in  KBr 


Size  range 
(pm) 


Peak  position. 


Mass  extinction  coeff, 
(em'/g) 


1000 

108G 

1084 

1087 


1,300 

8,000 

15,400 

33,100 


B.  Shape  effects 

Having  eliminated  the  size  dependence  problem  our  next 


The  implications  of  these  results  are  clear.  To  avoid 
the  messy  and  complex  problem  of  dealing  v;ith  extinction 
varying  in  a  complicated  v;ay  with  si.zc  distribution  ojie 
should  experimentally  separate  out  a  size  fraction  for  which 
this  problem  docs  not  exist.  Our  experiments  on  quartz 
particle  segregation  have  convinced  us  that  this  is  possible. 
The  technique  has  therefore  been  adopted  as  standard  practice. 
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concern  v/as  to  compare  measurements  with  the  theory  for  spheres. 
In  this  case  and  in  the  rest  of  our  v;ork,  absolute  normaliza¬ 
tion  of  the  extinction  coefficient  by  unit  volume  (rather  than 
by  unit  mass)  has  been  used.  It  is  important  for  us  to  point 
out  that  v/e  are  not  making  any  arbitrary  scaling  of  either  the 
experiment  or  the  theory,  as  is  ver_  frequently  done  in  such 
comparisons.  riguro  (7)  shows  the  volume- normalized  extinc¬ 
tion  coefficient  detcirmined  experimentally  from  the  small  size 
friiction  of  quartz  powder,  compared  v;ii.n  the  calculated  extinc¬ 
tion  based  on  accurately  measured  bul);  optical  constants  for 
quartz  (Spitzer  and  Kleinman  1961)  .  Tlic  agreement  is  seen  to 
be  very  poor.  Maximum  extinction  shov/s  a  disagreement  of  more 
than  a  factor  of  five,  and  widths  and  sliaptis  of  the  strong 
extinction  band  compare  poorly.  Tlie  probleni  is  evidently  shape 
effects.  VJe  have  used  a  theory  for  splieres  to  apply  to  highly 
irregular  particles  (see  inset  of  figure  4) .  Our  theory  block 
of  section  II  is  not  wrong,  but  it  is  inappropriate  to  the 
measurement  bloclc.  It  is  evident  that  use  of  tlie  sphere  theory 
to  invert  the  measurements  for  optical  constants  would  likely 
lead  to  very  bad  results.  One  needs  either  to  create  spheres 
of  the  particles  or  to  construct  a  tlicory  capable  of  a  better 
representation  of  irregular  particles.  Since  it  appears  im¬ 
practical  for  us  to  malce  spheres  (less  t]ian])jm  in  diameter)  of 
all  solids  v;e  v/ish  to  study,  our  alternative  has  been  to  attempt 
development  of  a  more  satisfactory  t)ieoretical  treatment. 
Fortunately,  for  particles  of  this  size  range,  tlie  Rayleigli 
theory  is  a  rather  good  approximation,  and  shape  effects  can  be 
easily  treated,  to  some  extent,  by  considering  ellipsoids  v;ith 
a  v/ide  variation  of  axial  ratios.  Details  of  this  development 
are  given  in  the  next  section.  In  the  lov;er  curves  of  figure 
(8)  V'/e  prov'iev/  our  theoretical  results  for  quartz  (one  of  tlie 
bettor  cases) .  Results  of  extinction  from  the  distribution 
of  Rayleigh  ellipsoids  calculation  arc  compared  with  the 
measured  extinction  curve,  with  no  arbitrciry  .scaling.  It  is 
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Figure  7:  Comparison  of  extinction  calculated  for  spherical 
particles  of  quartz  in  KBr  with  measured  extinction  for  a 
small  size  fraction  of  quartz  powder  in  KBr. 
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Figure  8:  Comparison  of  calculated  ani  measured  extinction 
for  quartz  particles.  The  upper  curves  compare  calculations 
for  a  continuous  shape  distribution  of  Rayleigh  ellipsoids  (CDE) 
with  calculations  for  spheres.  The  lower  curves  compare  the  CDE 
calculations  with  measurements. 
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obvious  that  the  accounting  for  of  shape  distribution  has 
markedly  improved  the  agreement  of  tlicory  and  experiment  and 
shows  promise  of  enabling  us  to  extract  optical  constants  for 
real  powders  composed  of  highly  irregular  particles. 

IV.  THEORY  OF  SHAPE  niSTRIHUTION  IN  THE 
IlAYLEICH  ELLIPSOID  APPROXJ  IVvTION 

The  average  absorption  cross  section  for  a  collection  of 
randomly  oriented  identical  homogeneous  ellipsoids  that  are 
sufficiently  small  compared  to  the  v;avclongth  of  the  incident 
light  may  be  v/ritten  in  the  form  (van  de  Hulst  19S7,  chap.  6) 

3  1 

<Cabs^  =  kv/3lm  {  I^(P  +  I-)"  ,  (8) 

6  -  (c  -  1)'^  , 

where  v  =  4iralDp/3  is  the  volume  of  an  ellipsoid  v;ith  serniaxes 
of  length  a>^b^c  ;  c  =  e,  tiCj  is  the  complex  dielectric 

function  of  the  particle  relative  to  that  of  the  surrounding 
medium,  and  k  =  2it/X  is  the  wave  number  of  the  surrounding 
medium.  The  L ^  '  s  are  geometriciil  factors  related  to  the  ratios 
of  semi-major  axes  for  the  ellipsoid  (van  de  Hulst  1957,  p.  71). 
Only  two  of  the  geometrical  factors,  which  we  shall  take  to 
be  L^  and  are  independent  because  of  the  relation 

t  Jj  ^  ^  ~~  i  * 

Suppose  now  that,  in  addition  to  being  randomly  oriented, 
the  collection  consists  of  all  possible  ellipsoidal  shapes; 
i.e.  the  georaetriccil  factors  are  not  restricted  to  a  single 
set  of  values  but  arc  distributed  according  to  some  shape 
probability  function  P(Lj^,Lp)  .  All  ellipsoidal  shapes  are 
represented  by  points  in  the  sluidcd  triangular  region  shown 


Figure  9:  Domain  of  definition  of  the  shape  probability 
function.  Outside  of  the  region  F  bounded  by  =  Of  L2 
+  ^2  ■=  1»  ^(Lj^,L2)  vanishes  identically. 
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'  in  figure  (9) .  However,  it  is  convenient  for  purposes  of 
integration  to  define  P  on  the  larger  triangular  region  T, 
which  is  composed  of  six  equal  area  regions,  each  of  which 
corresponds  to  one  of  the  six  possible  v/ays  of  choosing  the 
relative  lengths  of  the  ellipsoid  axes  a,  b,  c.  The  shape 
probability  function  is  normalized  to  unity  on  F ; 


II  dL2  =  1  (9) 

r 

The  absorption  cross  section  averaged  over  the  shape  distribu¬ 
tion  and  over  all  orientations  is,  therefore. 


<<C  .  >> 
abs 


=  ^(Lj  .bjUlLi  dLj  (10) 


=  kv/3  Im  ^3  ^ 


where 


= 


[  f 

J  J  3  +  L 


dLf  dL2 


=  11-^ 
B  + 


^(L;^,L2) 


dL2 


3  +  L, 


ff  ) 

- dL 

JJ3+l-bT-L„ 


dL, 


1  "2 


and  3  =  (c  -  1)~^.  The  integral  of  a  function  f(Lj^,L2)  over 
r  may  be  written  as  an  iterated  integral : 


II  f(Lj^,L2)  dLjL  dL2 
r 


1  1-Lj^ 

I  dbjL  I  f(Lj^,L2)dL2 
0  0 
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We  have  assumed  that  all  particles  have  the  same  volume  v; 
however,  if  there  is  no  correlation  botv/een  shape  and  volume, 
then  the  total  absorption  cross  section  of  the  collection  is 

7/1S5v>  +  '^2  '^3^ 

where  7/  is  the  total  number  of  particles  per  unit  volume  and 
<v>  is  the  average  particle  volume.  It  has  also  been  implicit¬ 
ly  assumed  that  ^(Lj^,L2)  is  continuous;  this  is  not  a  necessary 
restriction,  hov;evor,  and  wo  can  take  into  account  discrete 
distributions  by  replacing  the  above  integrals  vjith  summations 
over  the  discrete  set  of  points  in  1. 

Perhaps  the  simplest  conceivable  distribution  is  one  for 
which  all  shapes  are  equally  probable,  in  the  sense  that 
^(L^,L2)  =  2,  and  the  integrals  arc  readily  evaluated 

■^1  “  ^2  "^3 

Therefore,  the  average  cross  section  is 

>>  =  k<v>  Im  log  c)  .  (11) 

^  ^  c-1 

In  the  above  expression  log  denotes  the  pri ncipal  value  of 
the  logarithm  of  a  complex  number  z  =  rcl®  (Churchill  1960, 
p.  56)  : 

Other  functions  d’(L^,L2)  could  be  chosen  to  represent 
a  shape  distribution  of  particles.  The  one  chosen 
has  the  advantage  of  uniformity  in  the  bjbp  plane  of  figure  (9) 
and  simplicity  of  the  resulting  expression  (11) .  The  real 
test  of  its  practical  value,  of  course,  lies  in  how  well  it 
describes  th.e  properties  of  actual  particles. 
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V.  DIRECT  COMPARTEON  OF  EXPERIMENTAL  AND  THEORETICAL 
EXTINCTION  MEASUREMENTS 

Having  developed  a  theoretical  approach  for  treating 
nonspherical  particles  in  the  R£iyloigh  approximation  and  having 
successfully  developed  an  experimental  procedure  for  limiting 
the  particle  sixes  to  the  small  sixe  range,  wo  nov;  seek  to  find 
hov/  well  the  theory  and  experiments  agree  for  several  solids 
with  already  v;oll  knovm  optical  constants.  In  this  section, 
then,  v/e  compare  experimentally  determined  extinction  spectra 
v/ith  calculated  extinction  spectra  based  on  the  known  optical 
constiints.  This  is  the  direct  problem.  The  indirect  problem, 
which  is  our  ultimate  interest,  will  treat  the  optical  constants 
of  these  solids  as  unknovm  to  see  how  v/ell  our  inversion 
technique  can  determine  optical  constants.  This  will  be  the 
subject  of  section  VI. 

Four  solids  have  boon  selected  that  have  strong  absorption 
bands  in  the  8  to  20  ym  si:)Octral  r^ange.  The  main  criterion 
for  selection  of  those  "standard"  solids  is  that  their  bulk 
optical  constants  be  already  v;cll  known  from  conventional 
measurements  on  single  crystals.  The  solids  chosen  are  quartz 
(f;i02)  t  calcitc  (CaCO^)  ,  SiC,  and  NiO.  All  of  these  solids 
were  available  commercially  as  fine  powder.  NiO  is  the  only 
truly  isotropic  solid  and  was  chosen  partly  for  that  reason. 

Sic  is  iniaxial  but  only  slightly  so.  It  is  of  interest  because 
it  has  one  of  the  largest  calculated  extinction  peaks  near  10  ym 
of  any  solid  we  have  found.  Qvuirtz  and  calcitc  are  both  very 
common  minerals  that  are  anisotropic.  Calcitc  is  particularly 
anisotropic  and  was  chosen  primarily  for  this  reason. 

A.  Sample  preparation  and  measurements 

Samples  of  small  sized  powders  dispersed  in  KBr  pellets 
were  prepared  for  extinction  measurements  in  the  following  way: 

(1)  The  solid  v.’as  cither  acquired  as  a  fine  powder  or 
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bulk  material  was  ground  in  a  ball  mill  to  produce  a  fine 
powder.  Most  solids  are  available  from  the  manufacturer  as  a 
fine  pov;der  since  many  of  the  production  techniques/  such  as 
wet  chemical  precipitation,  give  rise  to  such  a  form.  All  the 
materials  reported  on  in  this  study  v;cre  in  fact  purchased  as 
powdered  material. 

(2)  The  pov;dcr  v;as  segregated  by  Stokes'  lav7  settling  in 
either  v;ater  or  alcohol.  Some  details  of  this  have  already 
been  given  in  section  III.  After  dispersal  of  the  powder  in  the 
liquid  the  dispersion  was  allowed  to  stand  overnight,  after 
which  the  top  portion  of  tlie  liquid  was  drav;n  off  and  dried  in 

a  glass  dish  on  a  hot  plate  for  a  sufficient  number  of  hours 
to  drive  off  the  v;ater. 

(3)  A  small  quantity  of  the  sample  powder  (of  the  order 
of  .5  milligram)  v/as  weighed  and  mixed  with  about  0.5  gm  of  KBr 
pov/der  by  agitation  of  the  mixture  in  a  small  glass  vial  with 

steel  balls  for  12  hours  or  more.  This  appears  to  be  necessary 
in  order  to  break  up  clumps  of  powder  and  to  insure  uniform 
dispersal . 

(4)  h  1/2"  diameter  pellet  of  KBr  and  sample  was  pressed 
from  the  pov;dcr  mixture  in  a  special  pellet  die  obtained  from 
Barnes  E}igincering  Company.  Ten  tons  of  force  were  used,  re¬ 
sulting  in  pellets  that  are  almost  transparent  in  visible 
light.  Transmission  spectra  (T  vs.  v;avclongth)  were  recorded 
with  a  Perkin  Elmer  model  137  spectrophotometer  with  a  pure 
KBr  pellet  in  the  reference  beo  .  to  correct  for  reflection  at 
the  KBr  surfaces. 

(5)  Volume  extinction  cocf  f  j  ciciits  (a)  wore  calculated 
from  the  expression 

T  =  exp  (~a^) 

where  a  is  the  mass  density  of  sample  per  unit  cross  sectional 
area  determined  from  the  weight  of  the  pov/der  and  the  dimension 
of  the  pellet.  The  bulk  density  of  the  solid  is  denoted  by  p. 


B.  Extinction  results  compared  with  theory 


Figures  (LO-13)  shov;  the  comparison  of  measured  extinction 
coefficients  v/ith  calculations  for  spheres  and  for  a  distribution 
of  Rayleigh  ellipsoids  using  equation  (11).  Optical  constants 
derived  from  bulk  solid  Tnoasurements  reported  in  the  literature 
v/ere  used  in  the  calculations.  References  to  the  published 
optical  constants  are  given  in  table  II. 

Table  II 

Sources  of  the  measured  optical 
constants  used 


Solid 


Reference 


quartz  (Si02) 
calcite  (CaCO^) 
a-SiC 
NiO 


Spitzor  and  Kleinmnn  (1961) 
Hellv/cge  et.  al.  (1970) 

Spitzer,  Kleinman  and  V7alsh  (1959) 
Gielisse  et.  al.  (1965) 


Generally,  the  distribution  of  ellipsoid  calculations  (CDE) 
come  much  closer  to  matching  the  measurements  in  magnitude 
and  shape  of  the  extinction  feature  tlian  do  the  sphere 
calculations.  In  the  case  of  quartz  the  agreement  in 
magnitude,  width  and  s]iape  of  the  bands  is  excellont.  For 
CaCO^  there  is  about  a  factor  of  two  discrepancy  remaining 
between  the  CDE  and  measurements,  but  the  improvement  over  the 
spherical  theory  is  substfintinl .  The  magnitude  agreement  be¬ 
tween  CDE  and  measurements  is  good  for  SiC  alt) ough  measurements 
show  a  rather  extended  tail  of  absorption  to  lower  frequencies. 
It  is  possible  that  this  tail  is  due  to  the  presence  of  other 
forms  of  Sic  in  the  powder  than  the  a-SiC  for  which  theoretical 
calculations  wore  made.  SiC  occurs  in  a  number  of  different 
crystalline  forms  (see  the  various  articles  in  O'Connor  and 
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Figure  10:  Comparison  of  measured  extinction  (dashed  curve) 
for  quarts  powder  in  KBr  with  calculations  for  spheres  and 
for  a  distribution  of  ellipsoidal  shapes  (CDE)  in  th 
Rayleigh  approximation. 


11:  Comparison  of  measured  extinction  (dashed  curve) 
ci'om  carbonate  powder  in  KBr  with  calculations  for 
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Figure  12:  Comparison  of  measured  extinction  (dashed  curve) 

for  Sic  powder  in  KBr  with  calculations  for  spheres  and 

for  a  distribution  of  ellipsoidal  shapes  (CDE)  in  the  Rayleigh 

approximation. 


CEXT/VOL  (CM-1)(X103  ) 

OOO  lO.OGO  20.aCD  30.000  40-000  50.000  60-000  70-000 


-25- 


,  Smiltons  1959) .  Our  information  supplied  by  the  manufacturer 
gave  no  idea  of  the  percentage  of  SiC  polytypes  present  in  the 
powder.  In  tlie  case  of  NjO  the  maximum  of  the  measured  extinc¬ 
tion  falls  about  half  v;ay  betv/een  the  sphere  calculations  and 
the  ellipsoid  distribution  calculations,  but  the  shape  and  width 
of  the  extinction  is  modeled  more  correctly  by  the  CDE .  The 
fact  that  the  peak  extinction  measurements  are  higher  than  tlie 
CDE  curves  in  the  cases  of  caJeite  and  NiO  probably  implies  that 
the  particles  arc  a  little  closer  to  being  spheres  than  to  the 
wide  shape  distribution  described  by  our  distribution  function 
of  equation  (  9 1 .  Indeed  it  is  reasonable  that  this  shape  dis¬ 
tribution  might  be  too  "drastic".  A  n\imerical  integration  over 
L  values  more  closely  surrounding  the  spliere  point  (Lj^=L2=l/3) 
in  figure  (9)  might  be  bettor.  Further  v/ork  along  those  lines 
is  planned  in  the  next  yecir's  work.  Because  of  the  intograbili- 
ty  of  the  shape  distribution  function  L^)  =  1  and  because 

it  does  scorn  to  give  a  rather  good  representation  of  extinction 
for  some  solids,  v;e  continue  to  uro  it  for  our  attempted 
inversions  of  the  data. 

VI.  INVERSION  OF  THE  EXTINCTION  DATA  TO  OIJTAIN 

OPTICAL  CONSTANTS 

The  proceeding  results  liave  shown  that  the  now  treatment 
of  particle  .shape  distribution  gives  substantially  better 
agreement  with  extinction  measurements  on  real  particles  than 
does  sphere  theory,  even  though  tlic  nev;  theoiy  is  not  entirely 
satisfactory  in  every  case.  We  nov;  use  the  measured  extinction 
spectra  for  our  model  solids  (the  daslied  curv’cs  of  figures  10- 
13)  as  input  for  calculations  to  determine  optical  constants. 

In  this  way  we  treat  the  optical  constants  as  unknov;n  for  these 
solids  and  evaluate  our  determinative  procedures  to  see  how 
accurately  the  derived  optical  constants  agree  v;ith  the  known 
optical  constants  for  those  materials. 
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'  A.  The  method 


The  process  uses  an  interactive  computer  with  which  the 
operator  con  repeatedly  adjust  the  input  parameters  into 
theoretical  calculations  of  extinction  until  the  calculated 
spectrum  converges  satisfactorily  close  to  the  measured  spec¬ 
trum.  The  computer  is  programmed  with  the  theoretical  equa¬ 
tions,  v;hich  are  combinations  of  the  Lorentz  oscillator  theory 
for  e(4)  follov/ed  by  the  extinction  theory  for  a  shape  dis- 
tirbution  of  Rayleigh  ellipsoids.  The  principal  equations 
are  the  Lorentz  multiple  oscillator  equation  for  complex  e 


c 


=  e  + 
00 
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0) 


2 

P»  j 


+  iYjW 


(4) 


and  the  shape  distribution  equation  (CDE) 


<<C  .  >>  =  k<v>  Im{  — log  c)  .  (11) 

abs 

Input  parameters  are  the  oscillator  parameters 

w  . ,  Cl)  .  ,  and  y  . 

3  P/D  D 

Although  it  might  at  first  seem  to  be  a  hopeless  task  to  var 
so  many  parameters,  a  little  experience  with  the  Lorentz  theory 
enables  one  to  make  exccllant  estimates  of  some  of  the  parameters 
in  the  beginning,  as  v;o  now  indicate.  First,  the  can  often 
be  found  from  handbook  or  literature  values  since  for  most 
insulating  solids  it  is  equal  to  the  sejuare  of  the  index  of 
refraction  in  the  visible  or  near  infrared.  This  quantity  can 
be  determined,  even  on  microscopic  particles  of  a  solid,  by 
varying  an  immersion  oil's  index  of  refraction  until  the  par¬ 
ticle  as  viev;ed  under  the  inicrc^scopc  ^ippears  to  vanish.  The 
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Oil  and  the  particle  then  have  identical  indices  of  refraction. 
A  discussion  of  this  technique  along  v;ith  an  extensive 
collection  of  n  values  is  found  in  Tlio  Particle  Atlas 
(McCrono  and  Delly  1973) . 

Next  one  must  determine  the  number  of  oscillators,  j, 
by  considering  the  number  of  resolvable  peaks  in  the  extinc¬ 
tion  curve  for  the  frequency  range  of  interest.  In  most 
cases  of  strong  infrared  extinction  there  will  be  only  one 
dominant  oscillator  near  the  region  of  interest.  This  is 
true  of  CaCO^  as  regards  tl)o  1500  cm  ^  region,  for  example. 

In  NiO,  one  major  oscillator  giving  an  extinction  peak  near 
480  cm”^  is  suggested,  although  further  v;enker  extinction 
appears  to  be  present  near  580  cm”^'.  In  the  case  of  quartz 
there  appear  to  be  two  main,  overlapping  extinction  bands, 
near  about  1090  and  1150  cm"'^.  These  examples  show  that  the 
first  estimate  of  j  is  not  difficult.  Estimation  of  w.  and  m 
is  greatly  aided  by  a  little  familiarity  witli  the  behaviour 
of  the  oscillator  equations.  Shape  broadened  extinction  bands 
for  tliO  extremely  intense  bands  of  concern  in  tljis  work  are 
well  confined  to  the  frequency  region  for  which  Cjis  negative 
(see  paragraph  4.3  of  the  review  article  by  Huffman  1977). 

The  low  frequency  limit  of  this  negative  Cj  region  coincides 
closely  v;ith  for  an  isolated  band.  This  enables  one  to 
make  a  good  estimate  of  Wj  imnediately .  The  high  frequency 
limit  of  the  extinction  band  occurs  at  the  point  where  Ej  =■  0. 
Since  absorption  is  lov;  in  this  region  one  can  reasonably 
neglect  the  damping  term  in  the  deno/ninator  of  equation  (4) 
and  set  e j  =  0.  These  approximations  are  indicated  below. 
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Solving  for  w  and  calling  it  a)(c=0) 

a)(c=0)  =  I  ^  <  +  ‘^0 

^«>  *  . 

The  only  oscillator  paramotors  loft  to  bo  chosen  are  the  damping 
factors  y.  Usually,  for  strong  bands  Yijj<<oi^  which  puts  some 
limits  on  y.  Beyond  this  it  is  probably  best  to  determine 
by  mating  a  few  trials  at  matching  the  measured  extinction  on 
the  computer.  Following  these  first  guesses,  any  of  the 
oscillator  parameters  may  need  to  be  corrected  and  other  less 
important  oscillators  added. 

A  brief  summary  of  the  techniques  suggested  for  making 
the  first  guesses  is  as  follows: 

(1)  Select  from  handbook  or  literature  values. 

(2)  Select  from  the  low  frequency  edge  of  the 

extinction  band. 

(3)  Estimate  tj  from  the  high  frequency  edge  of  the 

Ir 

extinction  band  using  equation (12) . 

(4)  Determine  tlie  remaining  parameter  y  by  computer 
iteration  to  match  the  data. 

Lest  the  reader  feel  that  this  is  a  terribly  subjective  and 
"hit  or  miss"  technique,  v;e  v;ant  to  point  out  that  this  is 
exactly  analogous  to  the  time-tested  oscillator  method  for 
determining  optical  constants  of  bulk  solids  from  normal 
incidence  reflectance  data.  Most  of  the  single  crystal  optical 
constants  rofcronced  in  table  II  have  been  derived  in  such  a 
way.  The  only  difference  is  that,  in  the  present  case,  the 
theory  incorporates  the  much  more  complicated  extinction 
for  a  distribution  of  shapes  rather  than  the  simple  and 
unquestionably  applicable  Fresnel  equation  (2)  used  in  the 
bvilk  studios. 
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'  Treatment  of  Anisotropy 

The  above  technique  has  assumed-  that  the  solid  is  isotropic, 
i.e.  that  its  optical  constants  are  the  same  for  all  polariza¬ 
tion  directions.  A  crystalline  solid,  hov/ever,  may  be  either 
isotropic  (all  directions  the  same) ,  uniaxial  (two  major  polari¬ 
zation  directions) ,  or  bi-axial  (three  major  polarization  direc¬ 
tions)  .  The  difference  these  cases  vjould  make  in  the  extinction 
expression  (11)  is  indicated  below: 

<<C^j^g>>/k<v>  =  Im  { )  isotropic  solids 

=  2/3  Im  {)  uniaxial  with  eJ^c 

=  1/3  Im  { }  uniaxial  with  E||  c 
or  bi-axial 

E|^c  and  e||s  denote  the  electric  field  orientation  perpendicular 
and  parallel  respectively  to  the  unique  axis  c.  For  purposes 
of  inversion  v/e  have  merely  assumed  that  the  optical  constants 
arc  isotropic.  This  could  result  in  an  error  as  large  as 
a  factor  of  throe  in  the  absorption  attributed  to  a  given  band. 

One  could  rather  assume  the  factor  2/3  so  that  the  maximum 
error  in  absorption  would  be  a  factor  of  2 .  It  should  be  pointed 
out  that  this  magnitude  of  error  does  not  necessarily  propagate 
into  the  optical  constants.  Further  work  should  be  done,  however, 
on  the  problem  of  how  to  select  the  factor. 

B.  Results  of  the  inversion  procedure 

Results  of  the  scheme  for  obtaining  optical  constants 
are  presented  in  figures  (14)  through  (17) .  Derived  dielectric 
functions  are  shown  as  dashed  lines  and  the  ^’true"  (measured 
for  single  crystals)  dielectric  functions  are  represented  by 
full  lines.  The  scales  of  the  figures  have  been  chosen  to 
emphasize  the  regions  of  most  consequence  in  determining  small 
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particle  extinction  properties,  although  this  emphasis  truncates 
the  extreme  excursions  of  Cj  and  In  the  case  of  quartz  the 

powder  results  compare  very  well  with  measured  bulk  constants. 
Crystalline  quartz  hcis  two  sots  of  optical  coiistants,  one  set 
for  the  ordinary  (O)  ray  and  one  set  for  the  extraordinary  (E) 
ray.  For  the  large  feature  near  1050  cm”^  the  measured  optical 
constants  are  very  similar  for  the  tv;o  directions  and  only  one 
set  is  shown  to  minimize  confusion  in  figure  (14).  Since  the 
optical  constants  in  the  tv/o  directions  are  so  similar  near 
*  1050  cm~^  the  crystal  acts  very  much  like  an  isotropic  crystal, 
which  our  fitting  procedure  has  assumed.  For  the  structure  near 
800  cm  there  is  more  difference  for  the  two  directions,  and 
the  pov.'dor  results  contain  structure  from  each  polarization 
direction.  The  silicon  carbide  results  are  very  close  to 
bulk  measured  properties  —  so  close  that  the  agreement  may 
contain  an  element  of  good  luck.  I]ov;ovcr..  the  true  optical 
constants  for  SiC  are  acurately  modeled  very  well  by  a  one 
oscillator  fit  with  very  lov/  y.  Since  our  procedure  has 
little  sensitivity  to  small  y'-»  only  the  two  parameters 
(Oj  and  need  be  cliosen  correctly;  those  luive  evidentiy  been 
chosen  quite  v.’ell  in  the  fitting  procedure.  The  NiO  results, 
figure  (17)  ,  are  not  so  favorahio  in  conipar.ison  with  measured 
values.  In  the  case  of  CaCO^  the  results  arc.  rather  good 
despite  tlie  fact  that  the  anisotropy  is  not  ttiJ;en  into  acci  unt  • 
in  the  inversion.  VJe  treated  the  situation  as  thougli  it 
v;as  an  unknown  for  vfnich  v/e  had  no  kno’wlcukje  of  whether  the 
band  near  1400  cm  is  present  in  one^  two  oi  thrc'e  of  the 
major  polarization  direcLions.  Tlx'.  bull;  optical  constants 
in  the  comparison  are  for  the.  0-ray  v.’h.ich  is  totally  responsible 
to  the  band.  Overall,  the  inversion  tochninue  for  dc!tormining 
optical  constants  of  pov.'ders  gives  results  which  compare  quite 
favorably  Vvfith  mcasin.-ed  bulk  values  for  the  cases  of  our 
four  comparison  solids. 
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.  VII.  COMMENTS  ON  THE  SEARCH  FOR  MAXIMUM  EXTINCTION 

Although  our  primary  goal  has  boon  the  development  of 
a  technique  for  determining  optical  constants,  the  process 
of  repeatedly  going  back  and  forth  between  calculations  and 
experiment  for  a  number  of  solids  has  provided  some  insights 
into  the  larger  problem  of  how  to  produce  maximum  extinction 
in  the  infrared.  Several  points  of  tliis  nature  v;hich  were  not 
fully  realized  by  us  at  the  beginning  of  this  work  will  now 
be  discussed  briefly. 


A.  Selection  of  the  resonant  frequency 


The  equation  (8)  for  extinction  by  a  Rayleigh  ellipsoid 
clearly  shows  that  a  resonance  of  the  sort  we  have  been  con¬ 
cerned  with  can  bo  made  to  occur  at  any  point  in  the  negative 
Cj  region  by  varying  the  relative  dimensions  of  the  ellipsoid 
to  change  the  L^.  The  resonant  condition  follov/s  directly  from 
minimizing  the  denominators  of  each  term  of  equation  (8) , 
giving  rise  to 


e,  =  -  1]  ^  (13) 

This  equation  is  also  discussed  in  some  detail  by  Huffman  (1977) 
To  exemplify • the  placement  of  the  ellipsoid  resonances  wc  have 
calculated  extinction  for  a  randomly  oriented  ellipsoid  of  SiC 
with  the  particular  l.^'s  of  .1,  .3,  and  .6.  The  results  are 
shown  in  figure  (18)  along  with  our  CUE  calculations  and  the 
real  part  of  the  dielectric  function  which  determines  the 
ellipsoid  resonance  position  according  to  the  above  formula. 

In  this  crystal  the  range  of  possible  resonance  frequencies  is 
from  about  780  cm”^  to  about  960  cm”^.  One  notices  that  the 
spread  of  the  shape  distribution  calculation  is  bctv;cen  these 
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rigura  18  <  Extinction  of  SiC  calculated  in  the  Rayleigh 
ellipsoid  approximation  for  randomly  oriented  identical  ellip¬ 
soids  with  shape  parameters  Lj  •  .1.  .3.  and  .6.  and  for  a 
continuous  distribution  of  ellipsoidal  shapes  (CDE) .  The 
real  part  of  the  complex  dielectric  function,  which  determines 
the  positions  of  the  peaks,  is  shown  below. 


limits  and  the  three  ellipsoid  peaks  fall  in  the  same  region. 

In  the  case  of  metallic  solids  it  is  possible  to  choose 
the  resonance  frequencies  over  a  much  wider  range.  The  negative 
Cj  region  for  many  metals  ranges  from  the  plasma  frequency, 
which  is  characteristically  in  the  far  ultraviolet,  through  the 
visible,  infrared,  and  microwave  regions  to  zero  frequency. 

Thus  a  resonance  in  a  metal  particle  could  conceivably  be 
designed  to  occur  at  any  frequency  in  the  infrared. 

A  calculated  example  of  extinction  for  one  particular  ellipsoid 
(Lj=  .01,  .3,  and  .69)  of  aluminum  is  shov/n  in  figure  (19). 

Once  again  the  CDE  calculation  and  the  curve  are  given  for 
comparison.  These  tv/o  illustrations  emphasize  the  important 
point  that  the  resonance  positionsof  a  single  ellipsoidal  shape 
can  bo  selected  if  one  can  control  the  shape  parameters;  how-^ 
ever,  a  distribution  of  shapes  inevitably  decreases  tlie  strength 
of  the  extinction  while  giving  a  broader  extinction  band.  It 
appears  to  be  technologically  impossible  to  produce  a  cloud  of 
sub-micron  particles  composed  of  identical  ellipsoids,  all  having 
exactly  the  same,  chosen  axial  ratios;  therefore,  it  is  perhaps 
more  important  to  consider  the  extreme  ellipsoidal  shapes  of 
spheres  (L^'s  =  1/3),  discs  1/2, 1/2,0)  and  needles  (Lj  = 

1,0,0).  These  shapes  can  frequently  be  taken  by  real  particles. 

B,  Practical  maximum  extinction  for 

spheres,  needles,  and  discs 

Since  the  heights  and  widths  of  the  extinction  peaks  are 
determined  by  optical  constants  of  the  various  solids,  not  many 
of  which  are  knov.Ti ,  it  is  not  possible  for  one  to  say  what  the 
maximum  possible  extinction  may  be  among  all  solids.  We  can, 
however,  choose  a  few  examples  from  promising  ionic  crystals 
having  measured  optical  constants  in  the  vicinity  of  their 
intense  absorption  bands  in  the  infrared.  Calculations  for 
peak  extinction  normalized  per  unit  volume  of  spheres,  needles. 
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Figure  19  J  KxtJnction  of  aluminum  calculated  in  the  Rayleigh 
cllipr.oid  approximation  for  rantlonly  oriented  identical  cllip- 
Boicls  with  ship*'  parainot  err.  a  ,01,  ,3,  .69  and  for  a 
continuous  disstribul  ion  of  eilipEoidal  shapes  (CM*),  The 
real  part  of  the  complex  dielectric  function,  v.*hich  determim'S 
the  i>o.Hitiono  of  the  p'^aKn,  in  rhown  Iflow, 
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'  and  discs  of  three  solids  are  listed  in  Table  III.  The  metal 
aluminum  is  included  in  the  three  for  comparison. 

The  first  important  point  to  observe  is  that  the  sphere 
wins  in  all  cases  in  the  competition  for  maximum  extinction 
among  the  three  limiting  shapes  in  the  vicinity  of  their 
shape-dependent  resonances.  The  magnitude  of  the  maximum 
extinction  for  spheres  is  determined  by  how  small  Ojis  at  the 
resonance  frequency.  For  these  examples  the  maximum  value 
of  extinction  is  about  4  x  10^  cm”^  in  the  vicinity  of  10  ym. 

For  the  present,  at  least,  this  represents  about  the  maximum 
expected  extinction  per  unit  volume  for  spheres  of  real 
material. 

C,  Practical  maximum  extinction  for  shape  distributions 

It  has  already  been  discussed  that  a  distribution  of  shapes 
gives  a  broader  but  less  intense  extinction  band  than  single 
shapes.  If  the  CDE  expression  is  applicable  to  real  particu¬ 
lates  one  can  determine  an  approximate  upper  limit  for  extinc¬ 
tion  by  a  collection  of  irregular  particles,  since  equation  (11) 
usually  tends  to  a  I’miting  value  in  the  vicinity  of  the  shape- 
dependent  resonances  with  low  damping.  The  reasoning  follows: 
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For  particles  in  air  the  limit  is  about  3  x  10  cm  for  a 
frequency  corresponding  to  about  10  ym  v/avelength.  This 
limiting  curve  is  shown  in  the  last  figure.  Our  highest 
extinctions  measured  in  this  v/ork  are.  consistent  with  this 


Average  Absorption  Cross  Sections  per  Unit  Particle  Volume  for  Randomly  Oriented 


limit  or  slightly  greater.  This  might  be  considered  to  be 
the  maximum  extinction  for  a  wide  distribution  of  shapes. 

The  considerably  larger  peak  extinction  value  of 
5  -1 

o  =  2.1  X  10  cm  has  been  reported  for  a  small  size  fraction 
of  powdered  talc  (Dorschner  et.al.l979),  a  mineral  with  the 
chemical  composition  Mg2  (OH)  2  (^^1205)  2  which  is  the  primary 
basis  for  the  common  talcum  powder.  Those  measurements  were 
discussed  briefly  by  Bohren  and  Huffman  (1979b)  and  are 
included  in  figure  (20).  it  appears  that  the  exceptionally 
large  extinction  values  may  be  due  to  the  particles  being 
nearly  spherical.  Since  talc  is  one  of  the  softest  of  all 
solids,  our  speculation  has  been  that  in  the  process  of 
forming  the  talc  into  a  KBr  pellet  the  high  pressure  molds 
the  talc  grains  into  an  approximately  spherical  shape,  con¬ 
siderably  enhancing  the  peak  extinction.  An  attempt  was  made 
to  use  our  technique  to  invert  the  extinction  data  for  talc 
using  the  ellipsoid  distribution  theory,  but  it  was  impossible 
to  decrease  y  enough  to  give  calculated  values  as  high  as  the 
measurements.  In  this  special  case,  because  of  the  softness  of 
talc,  it  proved  possible  to  press  a  pellet  of  pure  talc  with 
almost  bulk  density  for  specular  reflectance  studies.  Analysis 
of  the  specular  reflectance  curve,  treating  the  sample  in  this 
case  as  a  smooth,  bulk  sample  gave  optical  constants  for 

extinction  calculations.  The  results  were  that  the  CDE 
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calculations' peaked  at  about  2x10  cm  using  the  inferred 

optical  constants,  while  the  sphere  calculations  could  be 

5  • 

made  to  peak  at  about  the  measured  valueof  2.1x10  by  choosing 
a  value  of  the  damping  constant  within  the  uncertainty  of 
the  optical  constants  fit.  The  conclusion  from  this  approach 
agrees  with  our  earlier  conclusion;  that  is  that  talc  particles 
are  much  more  spherical  after  being  pressed  into  the  pellet 
than  most  collections  of  particles  in  KBr. 


C.  Comparison  of  some  extremes  in  infrared  extinction 

Our  last  figure  (20)  can  be  taken  as  a  sort  of  siimmary  of 
several  things  we  have  learned  in  relation  to  the  search  for 
maximum  extinction  in  the  10  pm  spectral  region.  A  number  of 
materials  and  limits  are  presented,  all  plotted  in  the  same 
voliime-normalized  extinction  units  with  no  adjusted  scaling. 

A  phosphorous  smoke  aerosol,  which  is  a  common  military 
screening  smoke,  has  been  plotted  from  the  published  measure¬ 
ments  of  Carlon  et.  al.  (1977) .  Calculations  on  SiC  for  spheres 
and  for  a  continuous  distribution  of  ellipsoid  shapes,  which 
have  been  shown  on  a  different  scale  in  figure  (12) ,  represent 
extinction  to  be  expected  for  spheres  and  for  a  wide  distri¬ 
bution  of  shapes.  Recall  that  measured  extinction  for  SiC 
particles  falls  near  the  CDE  calculations  although  the  band 
shape  is  somewhat  different.  Talc  measurements  are  included 
to  indicate  how  large  the  extinction  peak  may  be  in  some 
almost-spherical  silicates.  The  jiractical  maxima  for  spheres 
are  shown  as  a  series  of  dashed  peaks  to  remind  the  reader  that 
this  level  of  extinction  is  appropriate  only  for  narrow  band 
absorption  at  a  few,  but  not  all,  of  the  indicated  wavelengths. 
Placement  of  the  practical  maxima  for  spheres  is  based  entirely 
on  a  few  solids  for  which  we  know  the  resonance  position, 

and  for  talc.  Also  shown  in  figure  (20)  is  the  theoretical 
maximiara  for  a  distribution  of  shapes  based  on  the  saturation 
of  our  expression  for  a  continuous  distribution  of  Rayleigh 
ellipsoid  shapes.  A  curve  of  the  "practical  upper  limits" 
for  absorption  only,  deduced  by  Carlon  (1979) ,  is  included. 

It  does  not  include  effects  which  arise  from  strong  variations 
of  optical  constants  in  the  vicinity  of  an  absorption, 
which  has  been  the  main  emphasis  of  our  work. 

This  collage  of  extinction  information  is  meant  to  give 
an  overview  of  representative  maximum  values  that  might  occur 
in  the  10  ym  spectral  region.  Although  the  total  picture  may 
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Figure  20 s  Volume-normalized  extinction  for  a  collection  of 
flfevoral  highly  absorbing  materials  and  for  several  practical 
limits  discussed  in  the  text. 


still  bo  a  bit  confused,  one  fact  is  clear.  Spheres  of  certain 
materials  having  strong  vibrational  absorption  bands  in  this 
region  produce  very  intense,  but  quite  narrow,  extinction  bands. 
To  produce  broader  bands  it  appears  to  be  necessary  to  reduce 
the  maximum  extinction,  either  by  varying  the  shape  distribution 
or  by  using  solids  with  different  optical  constants. 

VIII.  CONCLUSIONS  AND  RECOI-IMENDATIONS  FOR 
FUTURE  V;ORK 

The  best  method  we  have  developed  for  deducing  optical 
constants  of  powdered  materials  is  to  segregate  a  sub-micron 
size  fraction,  embed  a  weighed  portion  of  the  powder  in  a  KBr 
pellet  for  transmission  measurements,  and  fit  the  measured 
results  with  a  theory  composed  of  Lorentz  oscillator  expressions 
and  Rayleigh  theory  for  extinction  by  a  distribution  of  ellip¬ 
soidal  shapes  programmed  into  an  interactive  computer.  Our 
conclusion  on  testing  the  method  for  a  series  of  solids  with 
known  optical  constants  is  that  the  method  is  a  qualified 
success  but  could  use  some  more  refining  and  evaluation. 

It  is  clear  that  shape  effects  in  the  region  of  strong  absorp¬ 
tion  are  extremely  important.  Since  every  powder  may  have  its 
own  unique  shape  distribution  it  is  difficult  to  construct  a 
general  theory  of  shape  effects  that  works  perfectly  for  all 
powders.  The  other  side  of  this  shape  distribution  sensitivity 
is  that  control  of  particle  shape  can  markedly  enhance  peak 
extinction  or  broaden  and  widen  the  band.  Understanding  of 
these  effects  are  important  in  the  ultimate  goal  of  being  able 
to  engineer  the  optical  properties  of  particulate  systems. 

In  future  work  it  is  suggested  that  more  of  the  same  kind 
of  close  comparison  between  optical  constants  determined  from 
single  crystals  and  from  powders  be  carried  out.  Although 
most  solids  cannot  be  obtained  as  single  crystals  there  are 
still  quite  a  few  silicate  minerals  that  can  be  acquired  as 


natural  single  crystals,  whose  optical  constants  have  not  been 
measured  in  the  infrared.  These  silicates  all  have  strong 
bands  in  the  important  10  ym  region  because  of  their  common 
Si-0  vibrational  bands.  Measurement  and  cataloguing  of  such 
single  crystal  optical  constants  in  the  5-15  ym  region  coupled 
with  simultaneous  efforts  to  extract  optical  constants  from 
powders  of  the  same  materials  will  be  pursued  with  support  of 
next  year's  NASC  contract. 
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